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Abstract 

We have conducted a careful study of current- voltage {I-V) characteristics in fully integrated 
commercial PbZro.4Tio.6O3 thin film capacitors with Pt bottom and Ir/Ir02 top electrodes. Highly 
reproducible steady state I-V were obtained at various temperatures over two decades in volt- 
age from current-time data and analyzed in terms of several common transport models including 
space charge limited conduction, Schottky thermionic emission under full and partial depletion 
and Poole-Frenkel conduction, showing that the later is the most plausible leakage mechanism in 
these high quality films. In addition, ferroelectric hysteresis loops and capacitance-voltage data 
were obtained over a large range of temperatures and discussed in terms of a modified Landau- 
Ginzburg-Devonshire theory accounting for space charge effects. 
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I. INTRODUCTION 



It has been known for a long time that the description of ferroelectric oxides as insu- 
lators becomes less and less valid as their dimensions are reduced. The conductivity of 



thin films is often modelled usiner well estab 
„ew devices such as .es,stiven,en».esnB 



ished semiconductor theories and concepts for 
are already emerging which exploit the semi- 



conducting nature of thin ferroelectrics. The frequently observed deformation of hysteresis 
loops by leakage currents, which masks the true ferroelectric properties of the system, apart 
from being a nuissance to researchers, is not directly a problem for non-volatile memory 
applications. However, charge injection associated with finite conductivity is believed to 
be involved in accelerating fatigue as well as raising other reliability issues. The increased 
power consumption and heating due to leakage are also undesirable. Leakage measurements 
are therefore a crucial part of any electrical characterization of a ferroelectric device. These 
are often performed using automated systems designed for semiconductor characterization. 
However, as shown in this article, great care should be taken if one wants to obtain true 
steady state current- voltage (I-V) data due to the the much slower dielectric relaxation 
in ferroelectrics compared with usual semiconductors as well as resistance degradation at 
high fields. Such artefacts in I-V measurements can lead to incorrect interpretation of the 
conduction mechanism. 

This paper reports a detailed study of the /- V characteristics of Samsung's lead zirconium 
titanate (PZT) capacitors. Together with capacitance- voltage (C-V) measurements and 
ferroelectric hystesis (P-V) loops, the data are analyzed in the light of various conduction 
models with the aim of understanding the leakage mechanism. 



II. EXPERIMENTAL 



The 144nm PbZro.4Tio.6O3 films were deposited by a sol-gel method at Samsung, using a 
high purity (99.9999%) precursor, on a [111] oriented Ft bottom electrode, resulting, after 
annealing, in an almost epitaxial interface between the metal and the columnar PZT grains. 
For each capacitor the top Ir/Ir02 electrode was patterned as a two-dimensional array of 
interconnected ~13xl3yum^ squares with a total area A = 5.31 x 10~^cm^. The capacitor 
structures were fully integrated and embedded in Si, with various metallic layers connecting 
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FIG. 1: Typical current transients obtained at 42° C at various biases. Note the change of scale 
from logarithmic at low field to linear at higher biases. 

the top and bottom electrodes to the external tungsten contact pads. Current-voltage-time 
(I-V-t) measurements were performed using the Agilent 4155C semiconductor parameter 
analyser, dynamic capacitance measurements were obtained using the HP4192A impedance 
analyzer and the Radiant Technologies Precision Pro tester was used for ferroelectric hys- 
teresis loops. 



III. RESULTS 



A. Current-voltage characteristics 

The I-t response at 42°C of our capacitors upon a stepwise application of a constant 
bias is shown in figure ^ and is typical of leakage behaviour of ferroelectric thin films. At 
short times the behaviour is dominated by dielectric relaxation which often follows the 
empirical Curie- von Schweidler law (/ oc with < n < 1). This is then followed by a 
region of steady state leakage (in this case most clearly visible for the intermediate bias), 
which may or may not be observed depending on the point at which the current starts to 
increase again due to time- dependent dielectric breakdown or resistance degradation. It is 
common practice when measuring /- V characteristics to choose a waiting time (between the 
application of the bias and the current measurement) which is longer than the dielectric 
relaxation time at intermediate biases. However, as is clear from figure d this can lead to 
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an overestimate of the current at high biases due to resistance degradation and at low biases 
due to dielectric relaxation, sometimes leading to apparent negative differential resistivity 
or ohmic-like behaviour. 

To obtain true steady state current readings and avoid unnecessary stress to the sample, 
we performed I-t measurements at each voltage until a clear plateau in the current was 
reached. In a few cases at high biases and high temperatures where no clear plateau could 
be resolved before the onset of resistance degradation, the minimum value of the current was 
used. The tungsten pads were contacted using lO/im diameter tungsten probe tips, which 
were found to be more effective in achieving a good contact and penetrating through the 
thin tungsten oxide surface layer than thinner tips. Due to the high homogeneity of the 
films and the almost perfect reproducibility of all electrical measurements, the quality of 
contact could be checked by looking at the hysteresis loops, and in particular at the coercive 
voltages, which were found to be very sensitive to any contamination of the tungsten pad 
surfaces. Prior to taking any measurements the samples were poled at three times the 
coercive field (for ~60s) and at the highest measurement bias (for 10-20s, depending on 
the temperature) to avoid any diplacement currents and internal field changes due to the 
reversal of ferroelectric polarization when measuring the I-V characteristics. In this way 
the forward and reverse bias /- V measurements were obtained for a range of temperatures; 
we will define forward bias as the case when the top Ir02 electrode is at a higher potential 
than the bottom Pt electrode. 




Bias (V) 

FIG. 2: Forward (open circles) and reverse (filled circles) bias I-V at different temperatures. 
Figure 121 shows that at intermediate temperatures the /- V are quite symmetric apart from 
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the very low field region below ~ IV. At low and high temperatures, however, some asym- 
metry is seen. Although different capacitors were used to obtain the forward and reverse 
bias data, this cannot account for the observed asymmetry due to the high reproducibility 
of the measurements over many samples, especially in forward bias. The asymmetry is quite 
low, well below an order of magnitude for almost all biases. It is also interesting to note 
that whereas in low temperature measurements the reverse bias currents are higher, they 
are lower at high temperatures; this is not expected if the asymmetry is due to differences 
in the interfaces, such as different barrier heights, at the top and bottom electrodes. In 
general, reverse bias currents were somewhat less reproducible and less stable, with break- 
down occuring at lower biases. For this reason we focused more on forward bias data for 
the quantitative analyses, althogh some reverse bias I-V will also be discussed. 

In the hope of understanding the mechanism responsible for conduction in these films 
we have looked at a number of models commonly used to describe carrier transport in 
ferroelectrics and tried to fit their predictions as quantitatively as possible to our own data. 
In the following subsections we therefore present the same set of results analyzed in terms 
of these different models. 



Space charge limited (SCL) conduction 

When the internal field is dominated by space charge (either from free or trapped carriers) 
the current density J is, in general, expected to have a power-law dependence on voltage, 
J (xV^. In the case of discrete traps a quadratic voltage and cubic thickness {d) dependence 



is expected 
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where is the carrier mobility, e the relative dielectric constant and 6f the fraction of injected 
carriers that are free (i.e. not trapped); in the trap-free limit 6f = 1. Unfortunately, since 
these are commercial Samsung films of a fixed thickness, we have to rely purely on the 
voltage dependence when analyzing our data. 

From figure El it is possible to identify, for reverse bias data in particular, regions which 
look linear on the log-log plot and could therefore satisfy the power-law dependence. At low 
negative biases the room temperature exponent n was found to be 2.1, which is compatible 
with the trap-free or discrete trap case. However, an unexpected temperature dependence 
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was found, with n decreasing steadily to 1.65 at the highest temperature. The high-field 
exponents were in the range 3.0-4.8 for reverse bias and 3.5-6.9 for forward bias data respec- 
tively. Such high exponents are not common to SCL currents, and are certainly not possible 
for discrete traps. Values of > 2 are possible if traps are distributed within the bandgap 
and they are expected to be temperature dependent [3|. In the case of an exponential trap 
distribution characterised by a temperature Ti 

n — 1 = — 
T 

and for a Gaussian trap distibution with a standard deviation at we expect 

X o 277(7? 

(n- 1)2 = 1 ' * 



16PT2' 

However, neither of these relationships were well satisfied by our data. We therefore believe 
that SCL conduction is not likely to be the dominant leakage mechanism in our films, though 
other authors ^ ^ have tried to explain their data using this model. 



Schottky thermionic emission 

The slight asymmetry between forward and reverse bias I-V suggests that the interface 
properties might play a role and we therefore consider Schottky thermionic emission as a 
possible conduction mechanism. In this case the sample is considered as two back-to-back 
Schottky diodes, with the current controlled by the field-induced lowering of the orginal 
zero-field potential barrier height (pBo at the reverse (or forward) biased electrode when 
electrons (or holes) are the majority carriers. The current density is given by j?! 

1 *n-\2 



J = A*T' exp{-q<f)Bo/kT) exp(/3 V^^) (1) 

with the Richardson constant A* = A7Te'm*k'^/h^ related to the effective carrier mass m*, 
and p = (e/A;r)(e/47reeo)^/^ k and h are the Boltzmann and Planck constants and q and e 
are the carrier and electron charges respectively. The important parameter in equation ^ is 
the field at the potential maximum which determines the barrier lowering. We consider 
two cases. 
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Case 1: uniform field (full depletion) 



Bulk ferroelectrics are good insulators and even in thin films it is common to assume 
that E = V/d. If the films are semiconducting this corresponds to the case of full depletion, 
which is expected for low concentrations of ionised dopants. Figure |21 shows that, especially 




FIG. 3: Forward bias data represented on the standard Schottky plot. The data which deviate 
fi'om the linear fit are shown in light grey. 



at higher fields, the data satisfy the Schottky equation quite weh; but we also have to see 
if the physical parameters that can be extracted from the data are sensible. The relative 
dielectric constant was calculated from slopes of the linear fits in figure IHl to be between 0.9 
and 1.1 over the temperature range shown. The physical range for e is limited from one side 
by the optical dielectric constant eoo ~ 6 J8| and from the other side by the static dielectric 
constant e^c, expected to be around 400 [2| for this composition along [111]. The calculated 
value of around e ~ 1 is clearly outside this range. From the temperature dependence of the 
intercepts the Richardson constant and effective carrier mass were found to be 2.4 Acm~^K~^ 
and 0.02me respectively, where mg is the free electron mass. These values also seem somewhat 
low compared to those previously reported for perovskites 3|- A barrier height of 0.9eV 
was also calculated, which is plausible from band diagram considerations. 



Case 2: -partial depletion 



In more semiconducting films, the dopant concentration N^) may be sufficiently high to 
allow the applied voltage to be dropped across a thin depletion layer near the reverse biased 



electrode (for an n-type or forward biased for a p-type material), resulting in a linear field 
profile in the sample and a modified value for Em which also depends on the built-in bias 

V edcCo V Q J 

The kT/q term originating from the free carrier contribution to the space charge within 
the depletion layer may usually be neglected. 

Recently Pintilie et al. fll^ have analysed their epitaxial single crystal PZT data using 
the partially depleted thermionic emission model, modified to account for the ferroelectric 
polarization, and have shown that they fit the log( J/T^) ~ + VuY^^ law. We performed 
a similar analysis on our films, fitting the I-V to the quarter-power law assuming various 
values of Vu and extracting the physical parameters from the slopes and intercepts of our 
fits. 

We found that the quality of the fit could be adjusted by varying V^, with the high 
temperature data fitting best for low or negative values of Vu and the low temperature data 
giving the best fits assuming a high positive built-in bias. From the temperature dependence 
of the intercepts of log( J/T^) vs. {V + VuY^^ plots, the effective carrier masses and barrier 
heights were extracted and are shown in table |l] for different assumed values of Vu- The 



Vu (V) 


m* /nie 


(eV) 


-0.05 


20 


1.34 





2.1 


1.28 


0.1 


0.10 


1.21 


0.2 


0.010 


1.16 


0.4 


0.00026 


1.09 


0.6 


1.27x10"^ 


1.02 



TABLE I: Effective carrier masses are highly dependent on the Vu value assumed for the fit. It 
should also be noted that, due to the exponential dependence of m* values on the intercept, they 
are very sentitive to small errors. 



effective carrier masses in perovskites are expected to be of order a few mg so the most 
physically plausible value for Vu is close to zero. The log J vs. {V + Vufl^ plot with Vu = 
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is shown in figure 0] and, it is clear that only the high temperature data gives a good fit in 
this case. 
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FIG. 4: Forward bias I-V fitted to Schottky thermionic emission model under partial depletion 
with zero Vbi- 

Provided the mobility of the carriers is not too low, so that they only spend a short time 
near the potential barrier maximum, e in equation ([H) should be replaced by eoo- Substituting 
expression (0) for Em into (H)) shows that the slopes in figure El allow the product Nn/e^^edc to 
be calculated. Assuming eoo = 6 and e^c ~ 400 gives a rough estimate for the concentration 
of ionised impurities Nd ~ 10^°-10^^cm~^. 



Poole-Frenkel conduction 



In insulating samples, field-enhanced emission of carriers from traps can lead to Poole- 
Frenkel currents, which have a functional field dependence similar to Schottky emission 
currents with 

J = aoEexp (^-^^ exp^PppVE), (3) 

where Et is the trap energy, Ppf = 2/? and cxo is the sample dependent zero-field conductivity. 
Poole-Frenkel conduction is a bulk mechanism and therefore it is usual to assume E = V/d. 
Thickness dependence can be very useful in discriminating between bulk and interface effects, 
however, as already mentioned such a study was not possible in our case. Thus to test the 
Poole-Frenkel model we have to rely purely on the fitting of the I-V curves. Figure El shows 
that the data fit equation Q extremely well over the whole two decades in voltage at almost 
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all temperatures. The relative dielectric constant extracted from the slopes was found to be 
6.3 at room temperature and in the range of 6.3-7.4 for the six curves shown with no specific 
temperature dependence. These values are very close to the optical dielectric constant. The 
trap energy was found to be l.OeV. 
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FIG. 5: Forward bias data showing a very good fit to the Poole- Frenkel model over two decades in 
bias. 




B. Capacitance- voltage 

The capacitance was measured as a function of dc bias using a lOOkHz 25mV tickle 
voltage at high and low temperatures with typical curves shown in figure IHl^a). The room 
temperature zero-bias dielectric constant extracted from the capacitance measurement via 
e = Cd/Aeo was found to be around 460. The a-axis and c-axis dielectric constants for bulk 
PbZro.4Tio.6O3 are 498 and 197 respectively j^. We beheve our samples to be preferentially 
[111] oriented and therefore the extracted dielectric constant is not far from the value of 398 
expected for perfectly [111] oriented bulk PZT. Figure IHl^b) shows a Curie- Weiss plot for the 
inverse zero-bias dielectric constant obtained from the capacitance measurements. From the 
gradient of the plot and equation a Curie constant C = 4.44 x 10^K~^ was calculated. 

If the samples were partially depleted, a change in bias would result in a change of the 
depletion widths and consequently of the depletion layer capacitance. The concentration 
profile of ionised impurities can therefore in principle be calculated from the C- V measure- 
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FIG. 6: (a) Capacitance- voltage curves at different temperatures, (b) Curie-Weiss plot for the 
inverse of the relative dielectric constant. 



ments using 



Nd{V) 



(4) 



where each value of V probes No at a particular depth within the film. In practice, however, 
only the mobile carriers respond to the ac tickle voltage used for the C-V measurements, 
thus only mobile carrier concentrations (n or p) can be obtained. Using this method and 
again assuming e^c = 400 we obtained a value of ~ 3 x lO^^cm^'^ for the mobile carrier 
concentration. 



C. Ferroelectric properties 

A typical ferroelectric hysteresis loop is shown in figure [7^, illustrating the very well 
defined coercive voltages and a symmetric shape with very little imprint. The remanent 
polarization of 23/iC/cm^ is significantly smaller than the 33yuC/cm^ expected for bulk [111] 
PZT of the same composition. The capacitors are known to be under compressive strain, 
so strain effects cannot be responsible for the lowering of the polarization as they would 
increase it. It was observed that Pf. progressively decreases after each etching step during 
the fabrication process and must therefore be associated with introduction of defects such 
as space charge or perhaps passive interfacial layers. 

Within the framework of Landau-Ginzburg- Devonshire (LGD) theory, for ferroelectrics 
with a second order phase transition and with polarization low enough that the terms 
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FIG. 7: (a) Room temperature P-V loop; (b) Square of the remanent polarization Pr as a function 
of temperature. 

in the free energy expression can be neglected, the hysteresis loops are given by a cubic 
equation 

E3 = aP3 + (3Pl (5) 

where the field and polarization are assumed to be along the c-axis (hence the subscript 
3)and a and (3 are the standard bulk LGD coefficients. Provided the coefficient coupling 
the polarizations along the different axes is smaller than f3 or at high enough fields, the field 
Elm] and polarization -P[iii] along [111] are also expected to be related by equation © but 



121 and 



the effect of substrate strain is 



with modified coefficients. As was shown by in 
also to rescale the LGD coefficients. From equation (0) the inverse susceptibility and the 
spontaneous polarization are given by 



a 



1 

eeo 



1 

X 



-2a 



{Tc - T) 
_ 2(Tc - T) 



(6) 
(7) 



The plot of P^ vs. T is shown in figure Efb) and was used to extract the values of (3' = 
1.52 X 10^ m^C~^F~^ and = 620 K; the latter together with the previously obtained 
Curie- Weiss constant C give the parameter a' as a function of temperature; the primes are 
used to show that these are no longer the bulk LGD coefficients, but their effective values 
modified by strain, s pac e charge and the fact that the samples are not [100] oriented. 



In another paper 



14| we have studied phenomenologically the effects of space charge on 



hysteresis loops for fully and partially depleted samples with various dopant concentrations. 
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Our results showed that increasing the amount of space charge caused a decrease in polar- 
ization and coercive field values as well as an increase in the observable dielectric constant. 
We also found that in the case of fully depleted films the effect of space charge was to 
simply rescale the original a and (5 coefficients; the P-E loops and the zero-field dielectric 
constant are completely depermined by these rescaled coefficients via (0) and ((Tj). This was 
not the case for partially depleted films, where the hysteresis loops no longer fitted a cubic 
relationship at all. 

The LGD fit using the experimentally extracted a' and coefficients is shown in fig- 
ure El^a) and is almost perfect, which would be consistent with the idea of a fully depleted 
sample, as would be expected for an insulating film. 



IV. DISCUSSION 



We have presented a large amount of data analyzed in various ways to try and find some 
clues as to the nature of the leakage mechanism in these commercial PZT capacitors. We 
now try to summarise our findings and make some conclusions. 

SCL conduction has been ruled out as the l-V do not seem to satisfy the required 
power law or the possible temperature dependences. We also rule out Fowler- Nor dheim 
tunnelling precisely because of the strong temperature dependence; the /- V also did not fit 
the log//l^^ oc XjV law for any significant voltage range. 

Schottky thermionic emission assuming a uniform internal field distribution yields an 
unphysical value for the dielectric const aiit, as does the modified Schottky model proposed 
by Simmons [15;] and used by Zafar et al. Ig, 12 1 to fit their data on (Ba,Sr)Ti03 thin films. 
Thus both these mechanisms can be ruled out. Fow completeness, we should note that 
the "standard" Schottky equation and that due to Simmons are just limiting cases of 
the more general thermionic emission- diffusion model by Crowell and Sze (see 0,0^) which 
predicts the same J-E expression as (0) but with A* replaced by A**{E). Due to the unknown 
field dependence of A** it is difficult to fit data to this generalized model, however, this field 
dependence is usually much weaker than exponential and thus any functional deviations 
from are significant mainly at low fields. Our data spans a large range of biases and we 
have already ruled out the two limiting cases of the thermionic emission-diffusion model, 
thus it is probably safe to rule it out completely for the case of uniform fields on the basis 
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of the unphysical value it gives for e upon fitting. 

If the sample is assumed to be partially depleted with the field at the potential max- 
imum given by (j2)) it is possible to obtain a fairly good fit to the Schottky equation (jT)) 
with physically plausible parameters, at least over a limited range of applied biases and 
temperatures. However, the fits yield a fairly high value for the concentration of ionised 
dopants that, according to our theoretical calculations 1^, is expected to produce a signif- 
icant deformation of the hysteresis loop, which is not observed, and is also in fact expected 
to modify the internal field profile from a linear to a cubic, making equation (0) invalid and 
resulting in the loss of self-consistency in the analysis. Even without resorting to theoretical 
modelling, there seems to be a large discrepancy between the ionised dopant concentration 
of 10^°-10^^cm~^ extracted from I-V data and the carrier concentration of lO^^cm"^ ob- 
tained from C-V measurements. The mobile carrier density which can follow the ac signal 
used in C-V measurements is expected to be lower than Nd ^ but not by three or four 
orders of magnitude. The partial depletion model, although unlikely in this case, is very dif- 
ficult to rule out definitively as it contains many adjustable parameters such as non- uniform 
dopant concentrations, unknown m*, Vu and the field dependence of A** which cannot be 
independently measured. 

The last possibility we have considered was Poole-Frenkel conduction which so far seems 
to fit all the experimental evidence available. The /- V give very convincing straight lines on 
the \og{I /V) vs V^^"^ plot for all temperatures and almost all two decades of applied bias. 
The extracted eoo ~ 6.3 shows good agreement with other independent measurements of 
the refractive index of PZT; and the assumption of a uniform field distribution is consistent 
with the shape of the P-E loops according to phenomenological modelling. The I-V also 
show little dependence on the polarity of the applied bias, again confirming that the leakage 
currents are probably bulk limited. 

Poole-Frenkel conduction in PZT has been previously reported by Chen et al. They 
observed highly symmetric /- V despite having different top and bottom electrodes (Au and 
Pt respectively) with Poole-Frenkel conduction identified above 40kV/cm by the extracted 
dielectric constant value of 6.3, which is in perfect agreement with ours. Nagaraj et al. 
[20I 1 also observed this mechanism for 120-360nm thin films with epitaxial (La,Sr)Co03 
electrodes. In their case, it was identified at fields above ~ lOOkV/cm and yielded an 
activation energy of 0.5-0.6eV attributed by the authors to ionisation of Ti'^+ to Ti^+. 
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Nagaraj et al. refer the 0.5eV activation energy to a private communication. Previous 



papers by the referenced author 
below the conduction band edge; t 
firmed in further pubhcations 



2l|, however, state that the Ti^+ centre hes at least leV 
le value of leV for the electronic trap depth was recon- 
23|. The ionisation energies of Pb ions 21| and oxygen 



vacancies 23[ are believed to be much lower compared with our value of Et leV, so 
electron trapping at Ti^+ sites may after all be responsible for the observed conductivity, 
as suggested in 20|. Et of leV is also similar to the activation energy for ionic transport 
of oxygen vacancies Q], however, the Poole- Frenkel model explicitly assumes a Coulombic 
attraction between a singly charged mobile and an oppositely charged stationary species; 
doubly charged carriers such as oxygen vacancies would change the slope of the log{I /V) vs 
V^^'^ plot. It is thus believed that we are observing electron or hole but not ionic conduction 
in these integrated PZT capacitors. 

Before concluding we would like to mention two more points. The first relates to dis- 
tinguishing Poole-Frenkel and Schottky conduction on the basis of the extracted dielectric 
constant. If the value extracted from the Poole-Frenkel plot agrees with the optical dielectric 
constant, while that from the Schottky plot does not, then it is likely that the mechamism 
is indeed Poole-Frenkel conduction. However, the reverse is not necessarily true. As was 

, in some cases a factor of k 

can enter the exponential term containing (3pp in equation giving equal slopes for both 
Schottky and Poole-Frenkel plots. 

Secondly, we note that our phenomenological arguments regarding the fitting of the P- V 
loop with coefficients obtained from measurements of e and Pr are only valid if the field is 
applied along the polar axis. If the field is along [111] (as in our case) and the polarization 
is along the c-axis, then the main contribution to the measured e will be from the larger 
a-axis dielectric constant e^, whereas the measured polarization will be the projection of P3 
onto the [111] direction. In this case the extracted Landau coefficients should not be able to 
reproduce the P-V loop. Thus is we assume that the fit in figure El^a) is not just fortuitous, 
then our polarization must be pointing along [111]. This may be possible, despite the 
tetragonal composition, through strain effects or through field induced change of symmetry 
to rhombohedral as described by Bell j26)]. However, so far we have not directly measured 
the direction of the polarization in our samples, nor have we explored in much detail the 
effects of different field orientations in our phenomenological modelling. Here we merely wish 
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to point out that P- V data may be of use in providing clues as to the extent of depletion in 
ferroelectric thin films. 

V. CONCLUSIONS 

The leakage mechanism in commercial PZT capacitors has been studied in detail. It 
was shown that standard automated /- V tests are not always appropriate for ferroelectrics 
due to the large difference in relaxation times for current transients at different fields. We 
have thus used current-time measurements to ensure that reliable steady state I-V were 
obtained. Together with C- V these data were carefully analysed and it was concluded that 
Poole-Frenkel emission of carriers from traps was the most likely conduction mechanism. 
A phenomenological model was used to confirm that the assumed field distribution was 
consistent with the observed shape of the ferroelectric hysteresis loop. 
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